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Abstract 

The analysis of nonlinear interaction of transversal electromagnetic field 
with Maxwellian collisionless classical and quntum plasmas is carried out. 
Formulas for calculation electric current in Maxwellian collisionless classical 
and quntum plasmas are deduced. It has appeared, that the nonlinearity 
account leads to occurrence of the longitudinal electric current directed along 
a wave vector. This second current is orthogonal to the known transversal 
current, received at the classical linear analysis. Graphic comparison of density 
of electric current for classical Maxwellian plasmas and Fermi—Dirac plasmas 
(plasmas with any degree of degeneration of electronic gas) is carried out. 
Graphic comparison of density of electric current for classical and quantum 
Maxwellian plasmas is carried out. Also comparison of dependence of density 
of electric current of quantum Maxwellian plasmas from dimensionless wave 
number at various values of dimensionless frequency of oscillations of electro¬ 
magnetic field is carried out. 

Key words: collisionless plasmas, Vlasov equation, Maxwellian plasma, 
Wigner integral, quantum distribution function, longitudinal electrical current. 

PACS numbers: 52.25.Dg Plasma kinetic equations, 52.25.-b Plasma pro¬ 
perties, 05.30 Fk Fermion systems and electron gas 

Introduction 

The nonlinear phenomena in plasma are stndied dnring the long time 
[I] - [I]- In work [Sj was the approach to stndying of dielectric fnnction 
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in quantum plasma is offered. Then in our works [9j-[[I2] was dielectric 
permeability in the quantum collisional to plasma is investigated. 

In works m and [m] plasma with any degeneration of electronic 
gas (plasma of Fermi—Dirac) was considered . It has been shown, that 
with use square-law decomposition of function of distribution is possible 
to reveal longitudinal electric current. This current is generated by the 
transversal electromagnetic held. Longitudinal current is perpendicular 
to the transversal current which is at the linear analysis. 

Let us notice, that existence of the longitudinal current generated by 
transversal electromagnetic held, it has been noticed more half century 
ago in work |2j. 

In the present work formulas for calculation electric current into 
Maxwellian collisionless plasma at any temperature (at any degrees of 
degeneration of the electronic gas) are deduced. 

It has appeared, that electric current expression consists of two sum¬ 
mands. The hrst summand, linear on vector potential, is known classical 
expression of electric current. This electric current is directed along 
vector potential of electromagnetic held. The second summand represents 
itself electric current, which is proportional to the square vector potential 
of electromagnetic held. The second current is perpendicular to the hrst 
and it is directed along the wave vector. Occurrence of the second current 
comes to light the spent account nonlinear character interactions of 
electromagnetic held with plasma. 

Then the kinetic equation with Wigner integral concerning quantum 
function of distribution is used. The nonlinear analysis is made and 
the formula for calculation longitudinal electric current is deduced. This 
current also is generated by transversal electromagnetic held. 

1. Vlasov kinetic equation and its solution 

Let us consider Vlasov equation describing behaviour of collisionless 
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plasmas 


df , df 


dt 


dr 


c 


df 


T^ + v^ + e E + -[v,H] 7 ^ = 0. 


dp 


(1.1) 


Vector potential we take as orthogonal to direction of a wave vector 


kA(r, t) = 0. 

in the form of the running harmonious wave 

A{r,t) = Aoe’**"-"*’. 


( 1 . 2 ) 


Scalar potential we will consider equal to zero. Electric and magnetic 
helds are connected with vector potential by equalities 

1 dA 

E =-—, H = rotA. (1.3) 

cdt ^ ^ 

The wave vector we direct along axis x: k = /c(l,0, 0), and vector 
potential of electromagnetic held we direct along axis y 

A = Ay{x,t){0, 1,0), Ay{x,t) - 


Then 

zc cJv 

Ay = - Ey, H = — Ey(0,0,1), 

(jJ (jJ 

cJ\ 

[V,H] = —Ey{Vy, -na;,0). 

CO 

Let us operate with method of consecutive approximations. Consi¬ 
dering, that the member with an electromagnetic held has an order, on 
unit smaller other members, let us rewrite the equation (1.1) in the form 


gf(k) Qj{k) 
-H Vx - 


+eEy 


' Qfik-d 

dpy 


1 - 


dt dx 

kv.r-\ kvydf^^~^^ 


CO 


+ 


CO dpj 


= 0, /c = l,2. 


(1.4) 


Here in zero approximation is the absolute Maxwell—Boltzmann 
distribution. 
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ks is the Boltzmann constant, T is the plasmas temperature. 

V p 

It is easy to see, that P = — = — is the dimensionless electron 

vt Pt 

velocity (or momentum), vt is the thermal velocity of electrons. 


Vt 


'2kBT 


m 


We notice that 


dm 


because 


5/(0) 


~ V. 


Therefore in hrst approximation Vlasov equation have the form 

a/(i) gfW _ Qfm 

+ Vx — 


dt dx ^ dpy 


(1.5) 


And in second approximation Vlasov equation have the following 
form 

5/(2) 5/(2) 


dt 


+ Va 


dx 


= —eE„ 


df^^^ u — kvx ^ kvydf^^^ 


dpy uj uj dpj 

We search solution in hrst approximation in the form 

In this approximation the equation (1.5) becomes simpler 

eEy dfo 


( 1 . 6 ) 


dfi ^ dfi 

From this equation we hnd 

-i{u - kvTPx)h = 


mvT dPy 


2eEy 

Pt 


MP)P, 


y> 
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from which 


fi = 


2ieEy PyfoiP) 


(1.7) 


Pt uj - kvrPx 

In the second approximation for function we search in the form 

/(2) = /(I) + /2 - El{x,t). 

Let us substitute in (1.6). We receive the equation 


df2 , df2 
- yVx - 




LU 


[UJ - kvE- -h kvy^ 

dpy dpx 


dt dx 

We transform this equation to dimensionless parameters 


df2 , j^df2 _ eEy 
+ VtPx-^ — 
ot ox UJPT 

From this equation we hud 


h P ^ I n P 

UJ - kvrPxj^r- + kvrPi 


dPy ' 


f2 = 


e^El 


ujpp[uj — kvrPx) 




dP„ 


h(P) 


dPx \uj — kvrPx 


( 1 . 8 ) 

We introduce new parameters q and Q, q is the dimensionless wave 
number, Q is the dimensionless frequency of elecromagnetic held, 


= A o = ^ 

^ kr ’ kxVT 


kT = ^ = 


Pt mvT 


h h 

Then we rewrite the equation (1.8) in the form 


/2 = 


e^E?. 


p‘^k‘^v^{fl — qPx) 


djPyMP)) I ,p2 d ( MP) \ 

dPy ^ '>dPPil-qPj 


. ( 1 . 8 ') 


Distribution function in square-law approximation on the held it is 
constructed 

/ = P’=/<“> + /l+/2. (1.9) 

where /i ,/2 are given accordingly by formulas (1.7) and (1.8). 


2. Density of electric current 
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Let us calculate current density 


'vfSv = euy / P/d^P = ev^ / P/id^P. 


( 2 . 1 ) 


By means of (1.8) it is visible, that the vector of current density has 
two non-zero components 

j {jxi jyi 0) • 

Here jy is the density of known transversal current, calculated as 


Jy = e4 / Pyf(fP = e4 / Pyhd^P 


( 2 . 2 ) 


This current is directed along electric held, its density is deduced by 
means of linear approximation of distribution function. 

Square-law on quantity of an electromagnetic held composed /2 the 
contribution to density of a current does not bring. Density of transversal 
current is calculated under the formula 


h ~ 




(2.3) 


mkr ^ — qP^ 

Here kr is the thermal wave number, 

Pt IJIVt 

^ ~ h ~ h ’ 

Let us calculate the longitudinal current. For density of longitudinal 
current according to dehnition it is had 

jx = e J v^fifv = e4 J Pxf2d^P- 

Having taken advantage (1.8'), from here we receive, that 


Jx 


k‘^nPQ. 


d{Pyfo{P)) , .2 A (J^\ 
dPy ^ ydP^\n-qP^J 


P^d^P 

D - qPx 


. (2.4) 


The hrst integral from (2.4) is equal to zero. Really, we will consider 
internal integral on Py 

\Py = + 00 


d 


dP 


iPyf0{P))dPy = PyfoiP) 


— 00 


= 0 . 


Py = — 00 
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Now in the second integral from (2.4) we will calculate internal integral 


on Pr 


d / MP) \ p,dP, 
dPx V n — qPx / n — qPx 


MP)P. 
(n - qP.W 


\Px=+00 OO 


m = -00 _r 


UP) p p^ 

Q — qPx — qPx 


= -n 


MP)dP, 
(Q - qPJ- 


Thus, equality (2.5) becomes simpler 


f MP)PyP 

J {Q-qP^Y' 


( 2 . 5 ) 


Double internal integral from (2.5) in plane (Py, P^) it is calculated 
in polar coordinates (Py = pcosc/?, P^ = p sin ip) 


2 ti oo 


e-Py-P" P^dPydP, = 


(y 2 O » « 

COS ipe~^ p difdp = —. 


—oo —oo 


Hence, the density of longitudinal current is equal to 


Jx 


e^E^Nq r e-^'dPx 

2 k‘^nPv^ypT J {qPx — D)^ 


( 2 . 6 ) 


Bbipancenne ( 2 . 6 ) npeflCTannivi b BHfle: 


4 °”® = Jci^.qWtrkEl. 


( 2 . 7 ) 


B (2.7) ai^tr ~ npoflOJiBHO-nonepennan npoBOflmviocTB, Jc(^5 q ) ~ ^es- 


pasMcpnan nacTb npoflojibnoro Tona, 


® r)2 7 ^ ^ eh / ujp 

prkr ^ p^\mv^) p^\kTVT^ ’ 



1 


00 


Jciyif q) — 




SnyPi J {qPx - 


Up - njiasMeHHasi (jiSHriviiopoBCKasi) HacTOTa, - SespasMepnasi njias- 
Mennasi nacTOTa, 


Up 




m 


n = 

a Lp 


k^vx 


Let us present the formula (2.6) in an invariant form. For this purpose 
we will introduce transversal electric held 


Etr = E 


kfEk 


E 


q(Eq) 


Now equality (2.7) we will present in coordinate-free form 
j‘“e = yn, qy,_t,.kEi = .U(l,qy,PlE,H]. 


The integral (2.6) is calculated with use of known rule of Landau as 
Cauchy type integral 


00 


'dr 


00 


= lim 


'dr 


{qr — e^o J i^qr — H + ieY 


-m- 


2q^ 


00 


+ V.p. 


^dr 


T=rt/q 


{qr — Uy 


On hg. 1-3 we will present graphics of behaviour of the real part of 
dimensionless quantity of density of electric current 


Re Jc(^) q) = V.p. 


1 


.-r- 


dr 


SiTyPr J {qr-Uy 


— 00 


On hg. 4-6 we will present graphics of behaviour of the imaginary part 
of dimensionless quantity of density of electric current, 
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Im Jc(Q, q) = 


.(e-Y' 




r=rt/q 






, n\2 


In case of small values of wave number from (2.6) it is received 


2a 


e^E^{x,t)N 

2m?‘uj'^ \(jj ) Sirmuj 


Q„\2 kEl 


nLp\^ i^xjy e nip 


n„\2 kEl 


'^Ltr 


kEl 


Q J Sirmcu Stt V Q 2 Qpxkr SttQ^ 

3. Quantum Maxwellian plasma 


( 2 . 8 ) 


Quantum Wigner distribution function was constructed in our work 
[7] in general case 

. , e«r^, /o(P + q/2)-/o(P-q/2) , eV^jFA)^ , 

^ ch uj — n^kP 2c^lE{(jj — n^kP) 

/o(P + q) -/o(P) /o(P-q)-/o(P) l /o(P + q) - /o(P - q) 

m — nTk(P + q/2) uj — nTk(P — q/2)J 4mc^/i to — n^kP 
By definition, the electric current density is equal 

j(r,Q = e J v(r,p,t)/(r,p,t)(i^n. 

Substituting in equality (4.1) obvious expression for speed 

eA(r, t) 


(3.1) 


v(r, P, t) = vtP — 


me 


and distribution function. 

Leaving linear and square-law expressions concerning vector potential 
of electromagnetic field, we receive 


J = e% 


vrPfi -A/o(P) 

me 




+ev 


T 


VTPf2 -A/i 

me 


Sp. 


(4.2) 
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Let us show, that the formula (4.2) for electric current density contains 
two non-zero components: j = {jx^jy^^)- One component jy is linear on 
potential of an electromagnetic held and are directed lengthways held. It 
is the known formula for electric current density, so-called "transversal 
current". The second component jx is quadratical on potential of held 
also it is directed along the wave vector. It is "longitudinal current". 

The hrst composed in (4.2) is linear on vector potential expression, 
and second is square-law. We will write out these composed in obvious 
kind 


; linear 



and 


. Vo(P + q/2)-/o(P-q/2) 
’ w - «rkP 


eA 

me 


MP) 




(4.3) 


:quadr 


= en 


T 


e^nTA(PA) 

mc^h 


[/o(P + f)-/o(P-p + 




e^4P(PA)- r /o(P + q) - /o(P) /o(P) - /o(P - g) 

Icj — nTk(P + q/2) cj — nTk(P — q/2) 


2c^ft^ 

e^nyPA" 

Amc^h 


[/o{P + q) - /o{P - q)] 


(PP 


uj — nykP 


(4.4) 


Expression (4.3) is linear expression of the electric current, found, in 
particular, in our previous work [7]. This vector expression contains only 
one the component, directed along the electromagnetic helds. Really, if 
wave vector to direct along an axis x i.e. to take k = A:(1,0,0), and 
potential electromagnetic helds to direct along an axis i.e. to take 
A(r,t) = {0, Ay{x,t),0), from the formula (4.3) we receive 


■linear 

Jy 


e^v^Ay 

mcq 


Here 


^ MP. + qm^MP. q/‘P) p2 ^ 

(4.5) 

’ kTVT 


fo{Px di q/2) 
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Let us consider expression for an electric current (4.4), proportional 
to square of potential of electromagnetic field. Let us notice, that the 
hrst composed in this expression is equal to zero. Hence, this expression 
becomes simpler 

iqua<ir = „,3 f\ eV^PiFAf p„(P + q) - /o(P) 

^ J 2c^h^iuj — utIcP) Vuj — n'rk(P + q/2) 

/o(P-q)-/o(P) l _ eVPAqo(P + q)-/o(P-q) 1 ^ 3 ^ 

w — nTk(P — q/2) J Amc^h cj — uykP 

Let us notice, that vector expression (4.6) contains one non-zero the 
electric current component, directed along the wave vector 

, <,uadr ^ f \\ MP. + q)-MP) . fo{P. - q) - MP) ] p 2 , 

^ 20^171^ J L qP^ + q‘^/2 — Q, qP^ — q^/2 — H J ^ 


/o(P, + <j)-/o(Px-g) 1 PJ^P 

2 qPx — n 


(4.7) 


Let us lead to the form convenient for calculations, the formula (4.7) 
for density of a longitudinal current. 

Let us consider the hrst integral from (4.7). We will calculate the 
internal integrals in a plane (Py, P^), passing to polar coordinates 

/„{P, ± q)P^dPydP, = „ 3 - , 

^ 2uyV7r 

Thus, size of generated longitudinal current into quantum plasma is 




equal 


4 f 

■quant _ _ V I 

^ AnPc^vry/^ J 

— 00 


_ e-T? ePPx-qf _ g-L? 

qP, + qy2 qP^- 12 -P 


+£-(^’.+9)’ _ e-(P.-9)" 


PxdPx 

qP^-Q' 


(4.8) 
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Let us transform the expression standing in integral (4.8). At first let 
us pass from potential to intensity of the held Ay = —(ic/uj)Ey. We will 
receive 

e^NAl _ e^NEl _ eulE^ 

16^^^/7TpT0J‘^ 

_ __ ^—E^ — ^ ^2 

167T^/npTkT^^q ^ 167T^/^TQ‘^q 

where the quantity of longitudinal-transversal conductivity ai^tr was 

introduced earlier: ai^tr = ^^‘1/{PThr) ■ 

Now equality (4.8) we will present in the form 

q)ai,tr‘^EyH,. (4.9) 

In (4.9) J{U,q) is the density of dimensionless longitudinal current, 

OO „ 

T (Q n)= ^ f - gV2]e-^ dr 

^ ’ IGtt^/ttQ J [{qr — Uy — q'^/4:][(qr — Uy — q'^] 

— OO 

At calculation singular integral from (4.10), which not writing out let 
us designate through /(n,g), it is necessary to take advantage known 
Landau rule. Then 



/(n, g) = Re /(n, q) ilm /(il, g). 


Here 


[3q^ + 2{qT — H)^ — g^/2]e 


Re/(n,g) = V.p. 


— OO 


[{qr — n)2 — ^^/4][(gr — n)^ — 
symbol V.p. means a principal value of integral. 


Im /(n, q) = —IT 


n 

q^ 1 q^ 


^-in/q+q/2p _ ^-iniq-ql2p 


+ 


+(w^ 


^-{n/q+qp _ ^-{n/q-qp 
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Equality (4.10) for density of longitudinal current we will present in 
the vector form 

jO"”* = J,(Q, q)<7i, trkEl, 


or 


Let us show, that at small values of wave number density longitudinal 
current both in quantum and in classical plasma coincide. 

According to (4.7) at small q after linearization /o(Pa;±g) and transition 
from vector potential to electromagnetic held we receive 


■quant _ 

Jx 


e^NqE 


eculkE^ 




ekEl /Q 


e^NqE, 


2m?VT^uj‘^ 


87rpTkT^uj‘^ STrpxkT^ V 

that in accuracy coincides with expression (2.8) for the classical plasma. 
These expressions we will copy in the vector form 

1 , _o 1 


kEl 


j'””' = -m.tr 




LO 


8ttQ^ c 

3. Conclusions 


E,H . 


In the present work the solution of Vlasov equation for collisionless 
Maxwellian plasmas is used. For the solution it the method of consecutive 
approximations is used. 

As small parametre the quantity of the vector potential of electromagnetic 
held (or to it proportional quantity of intensity of electric held) is considered. 

At use of approximation of the second order it appears, that the 
electromagnetic held generates an electric current directed along the 
wave vector, and proportional to the size of square of electric held. 
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Fig. 1. Real part of longitudinal electric current density, fl = 1. Curves 1,2 
correspond to classical Maxwellian plasma and classical Fermi—Dirac plasma, 
dimensionless chemical potential a = 0. 



Fig. 2. Imaginary part of longitudinal electric current density, D = 1. Curves 1,2 
correspond to classical Maxwellian plasma and classical Fermi—Dirac plasma, 
dimensionless chemical potential a = 0. 
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Fig. 3. Real part of longitudinal electric 

current density. Curves 1, 2 correspond to Maxwellian classical and quantum plasma. 



Fig. 4. Imaginary part of longitudinal electric current density. Curves 1, 2 
correspond to Maxwellian classical and quantum plasma. 
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Fig. 5. Real part of longitudinal electric current density of quantum Maxwellian 
plasma. Curves 1,2,3 correspond to values of dimensionless frequency of 
electromagnetic field 11 = 1.1,1, 0.9. 



Fig. 6. Imaginary part of longitudinal electric current density, of quantum 
Maxwellian plasma. Curves 1,2,3 correspond to values of dimensionless frequency 

of electromagnetic field 12 = 1.1,1, 0.9. 
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